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ABSTRACT 

The Crab Nebula has been observed by the HEGRA {High-Energy Gamma- 
Ray Astronomy) stereoscopic system of imaging air Cerenkov telescopes (lACTs) 
for a total of about 200 hrs during two observational campaigns: from Septem- 
ber 1997 to March 1998 and from August 1998 to April 1999. The recent detailed 
studies of system performance give an energy threshold and an energy resolution 
for 7-rays of 500 GeV and ~ 18%, respectively. The Crab energy spectrum was 
measured with the HEGRA lACT system in a very broad energy range up to 
20 TeV, using observations at zenith angles up to 65 degrees. The Crab data 
can be fitted in the energy range from 1 to 20 TeV by a simple power-law, which 
yields dVdE = (2.79 ± 0.02 ± 0.5) • 10-^(p|^)-2-59±o o3±o o5phm-2 s-^ TeV'^ 
The Crab Nebula energy spectrum, as measured with the HEGRA TACT system, 
agrees within 15% in the absolute scale and within 0.1 units in the power law 
index with the latest measurements by the Whipple, CANGAROO and CAT 
groups, consistent within the statistical and systematic errors quoted by the ex- 
periments. The pure power-law spectrum of TeV 7-rays from the Crab Nebula 
constrains the physics parameters of the nebula environment as well as the models 
of photon emission. 

Subject headings: TeV 7-rays: observations - Cerenkov telescopes: individual 
(Crab Nebula) -Supernova Remnant 

1. Introduction 

The Crab Nebula has been observed and studied over an enormously broad photon 
energy range embracing the radio, optical, and X-ray bands, as well as, high energy 
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7-ray region up to hundreds of TeV. The various theoretical scenarios of photon emission 
are primarily based on the Synchro Compton model (Gould 1965) which combines the 
synchrotron and inverse Compton (IC) emissions from high-energy electrons, which are 
accelerated up to ~ 100 TeV and interact with the magnetic field and the low frequency 
seed photons within the Nebula (de Jager & Harding 1992; Atoyan & Aharonian 1996; 
de Jager et al. 1996; Hillas et al. 1998). The predicted IC spectrum in the TeV energy 
domain appears to be very sensitive to the model parameters, such as the value of the 
magnetic field, the nature of the seed photons, the maximum energy of electrons, etc. 
Although IC scenarios of the photon emission are widely beheved to be appropriate for the 
Crab Nebula one can not exclude the possible contribution of 7-ray fluxes from 7r°-decay 
(see Atoyan & Aharonian 1996, Bednarek & Protheroe 1997). Thus, a measurement of 
TeV Crab Nebula spectrum sets major constraints on theoretical expectations and precise 
spectral measurements allow to fix the model parameters. 

The imaging air Cerenkov technique was successfully used for observations of the 
Crab Nebula in TeV 7-rays. Since the time of detection at a da confidence level by the 
Whipple group (Weekes et al. 1989) a number of observations of the Crab Nebula have 
been made at TeV energies (for reviews see Ong 1998; Catanese & Weekes 1999). By now 
the Crab Nebula is established as the standard candle of steady TeV 7-ray emission. The 
optical Nebula of the Crab has an angular extension of about 6 arc min (Hester et al. 1995). 
The standard analysis for the HEGRA system of lACTs provides an angular resolution 
of about 0.1° (see Figure 1). In the present analysis the Crab Nebula was assumed to be 
as a 7-ray point source (sec Figure 2). The detailed mapping of the TeV 7-rays from the 
Crab Nebula undertaken by Hofmann (1999) provided an upper limit for the Crab angular 
extension in TeV 7-rays of ~ 1.5 arcmin. No pulsed emission has been seen so far from the 
Crab (Gillanders et al. 1997; Burdett et al. 1999; Aharonian et al. 1999a). Vacanti et al 
(1991) provided a first measurement of the energy spectrum of TeV 7-ray emission from the 
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Crab Nebula using 10 m Whipple telescope data and derived a power-law spectrum index 
of 2.4 and a differential 7-ray flux at 400 GeV of 2.5 ■ 10~^° photons cm~^ s~^. Note that the 
absolute calibration of individual imaging air Cerenkov telescopes is quite difficult because 
there is no "test" beam of TeV 7-rays. The uncertainties in absolute calibration propagate 
to uncertainties of estimated 7-ray fluxes. Different methods of telescope calibration have 
been developed recently in order to reduce systematic errors which could influence the 
estimate of the absolute 7-ray flux and the slope of energy spectrum (see e.g., Prass et 
al. 1997). The energy spectrum measurements heavily rely on Monte Carlo simulations of 
the telescope response. Thus in the past the analysis of the same observational data using 
various methods of telescope calibration revealed, not infrequently, very different estimates 
of the telescope energy threshold and of the TeV 7-ray fluxes. These estimates may vary by 
more than a factor of 2. Thus, based on data taken with the prototype HEGRA imaging 
Cerenkov telescope (CT2), Konopelko et al. (1996) detected the TeV 7-ray signal from 
the Crab Nebula with a signal-to- noise ratio of lOcr and derived from the data a power- law 
index of the energy spectrum of 2.7 ±0.1 and an integral flux of TeV 7-rays above 1 TeV 
of 8(±l)stat(±2.4)syst • 10~^^ photons cm~^s~^. However, later on, the detailed treatment of 
the telescope hardware (optical smearing, photon-to-photoelectrons conversion efficiency, 
etc) allowed a more precise estimate of the TeV 7-ray flux using the single telescope data, 
as 1.5(±0.2)stat(+l-0 - 0.5)syst • 10"^^ photons cm-^s"^ above 1 TeV (Petry et al. 1996). 

Significant improvements in the telescope hardware as well as in the simulations and 
the data analysis (Fegan, 1997) recently provided measurements of the Crab Nebula energy 
spectrum over the energy range from 200 GeV up to 50 TcV by several groups, Whipple 
(Hillas et al. 1998), CANGAROO (Tanimori et al. 1998), and CAT (see Barrau 1998) 
exploiting the imaging Cerenkov technique. An additional measurement was provided by a 
low energy air shower array (Amenomori et al. 1999). However the uncertainties of the 7-ray 
flux estimates as well as of the energy spectrum slope remain rather large, in particular at 
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the high energy end of the measurements i.e. beyond 10 TeV. For example the Tibet data 
(Amenomori et al. 1999) show significantly higher (by a factor of 2) 7-ray fluxes compared 
with the results obtained using the lACTs in the energy range from 3 to ~ 18 TeV, and 
are in favor of a gradual steepening of the spectral slope at high energies. Thus additional 
precise measurements of the Crab Nebula energy spectrum in TeV 7-rays are of great 
importance. The HEGRA system of lACT provides such data. It was primarily designed 
for detailed spectral measurements in the TeV energy domain utilizing the advantages of 
stereoscopic observations. Stereo imaging gives several advantages for spectral studies, 
compared to a single telescope: (i) direct measurement of the shower impact parameter with 
an accuracy better than 10 m (ii) good energy resolution of 18% (in) wide dynamic range 
from 500 GeV to 20 TeV (iv) extended abilities for systematic studies using several images 
for an individual shower. The detailed systematic studies for the spectrum evaluation 
technique have been recently made using Mrk 501 1997 observational data (Aharonian et al. 
1999b, 1999c). The performance of the system was discussed by Konopelko et al. (1999a). 
Here we present Crab Nebula data taken with the HEGRA telescope system in 1997/1998 
and 1998/1999 observational campaigns. The data were analyzed using a new technique of 
energy spectrum evaluation for the stereoscopic observations. We also discuss the physics 
implications of the present results for the modeling of the TeV 7-ray emission from the 
Crab Nebula. 

2. Observational data 

The Crab Nebula was extensively observed with the HEGRA lACT system in two 
observing seasons from September 1997 to March 1998 and from October 1998 to April 1999. 
The observations were made with the stereoscopic system of lACTs which are located on 
La Palma, Canary Islands (Aharonian et al. 1999c). Each of the telescopes consists of a 
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8.5 reflector focussing the Cerenkov light onto a photomultipher tube camera. The 271 
photomultiphers in the camera were arranged in a hexagonal matrix covering a field of view 
with a radius of 2.15°. The telescope camera was triggered when the signal in two next 
neighbors of the 271 photo multiplier tubes exceeded a threshold of 8 photoelectrons, and 
the system readout started when at least two telescopes were triggered by Cerenkov light 
from an air shower. The detection rate was 12.6 Hz near the Zenith in December 1997 and 
dropped down to about 10 Hz in December 1999 for the 4-telescope system due to aging of 
the PMTs and reduced mirror reflectivity. 

The Crab Nebula was observed in a "wobble mode"; i.e., the telescopes were pointed in 
Dechnation ±0.5° away from the nominal Crab Nebula position (the sign of the angular shift 
was altered from one run of 20 min to the next). This is useful for continuous monitoring 
of the cosmic-ray background because it positions the OFF-source region symmetric to the 
camera center, and 1° apart from the ON-source region. Observations of the Crab Nebula 
at zenith angles up to 50 degree were made with 4 telescopes from 1997 September 1 to 
1998 March 29, for a total of 82.5 hr of data taken at good weather. Through the fire at 
the HEGRA site one of the telescopes was damaged and was out of operation for a month 
in October-November 1997. At that time. Crab Nebula observations were made with only 
three telescopes in the system, providing an event rate of 10.3 Hz near the zenith. Due to 
unstable weather and a substantial amount of dust which came from Sahara desert to the 
island, the average detection rate in 1998 February-March was reduced down to 10.7 Hz 
in observations near the zenith. At the beginning of 1998/1999 observational period 4 
telescopes were operational. Since October 1998 the HEGRA collaboration operates 5 
telescopes. However, for technical reasons one telescope (CT2) was out of operation since 
December 1998 until the end of the first observation campaign (April 1999). During the 
last observational period (since August 1999) observations at zenith angles less than 50° 
were taken for about 76.1 hrs. In addition, the observations at large zenith angles (LZA) 
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(50° < 6 < 65°), were carried out for a total of 24 hrs in order to study the performance of 
the telescope system at LZA and to extend the measurements of the Crab Nebula energy 
spectrum beyond 10 TeV. The total exposure times for the three periods are summarized 
in Table 1. 

Only data taken under good weather conditions were used in the analysis. In order 
to exclude data taken under less than optimal telescope performance conditions the entire 
database has been checked very carefully as follows: First, each night the compressed 
protocols of the system performance were transfered to one of the collaboration host 
institutes where they were scanned by software tools which closely monitor the status of 
the telescopes' hardware (single pixel rates, trigger rates of system telescopes, tracking 
accuracy, etc). This information was accumulated in a corresponding database which was 
used afterwards for a standard data reduction procedure. The final condensed data file for 
each particular run contains all the information needed for data analysis. In addition a 
specific software tool was developed which allows to control a posteriori for each data run 
(i) the system trigger rate, taking into account the zenith angle dependence (ii) the angular 
shape of the cosmic ray images, tested by a x^-criterion for the deviation of the mean scaled 
Width distribution for a single run from the corresponding average distribution filled over 
an extended sample of runs (Hi) the flatness of the ^^-distribution for the isotropic cosmic 
ray images over the full field of view (iv) the image Size distributions for each individual 
telescope. 

3. Analysis 

The stereoscopic imaging analysis of the data is based on the geometrical reconstruction 
of the shower arrival direction and the shower core position in the observation plane, as 
well as on the joint parameterization of the shape of the Cerenkov light images. The 
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simultaneous registration of several (>2) Cerenkov light images from an air shower provides 
an angular resolution of ~ 0.°1 for 7-ray showers. For each individual shower, stereoscopic 
observations allow to determine the position of the shower axis. Thus, at first only air 
showers within a certain impact distance Rq from the center of the telescope system were 
selected. The limiting upper radius of Rq = 200 m was used for zenith angles less than 50 
degrees, and a significantly larger radius of 400 m for the large zenith angle observations 
(> 50 degrees). The effective collection area in observations at LZA dramatically increases 
at high energies, far beyond the hmiting radius of 200 m (Konopelko et al. 1999b). For 
the data taken at zenith angles up to 50° an orientation cut < 0.05 [deg^] was applied, 
where 6^ is the squared angular distance of the reconstructed source position from the true 
source position. In addition the data were analyzed using the mean scaled Width parameter, 
< w >. To compensate for the dependence of the image shape on primary shower energy 
and distance from shower core to the telescope (impact parameter), the standard parameter 
Width (Fegan, 1997) (w*^), calculated for each telescope, is scaled according to the Monte 
Carlo predicted values, for 7-rays < w >^j, taken for the corresponding bin of reconstructed 
distance from the telescope to the shower core (i) and for the corresponding bin of image 
size (total number of photoelectrons in the image) (j) (Aharonian et al. 1999b, 1999c). The 
mean scaled Width parameter is defined for each individual shower as follows 

N 

<w>= 1/N ^ wV < w >J,. (1) 

k=l 

where is the number of triggered telescopes. This parameter was introduced in order 
to provide an almost constant 7-ray acceptance over the dynamic energy range of the 
telescope system. The optimum cut on mean scaled Width is about 1.1, which gives a 7-ray 
acceptance of ~ 60% at Small Zenith Angle (SZA). However, for a precise determination of 
7-ray spectra, a loose cut on mean scaled Width (< w >< 1.2) has been so far used in the 
data analysis in order to maximize the 7-ray acceptance and to minimize systematic errors 
related to cut efficiencies. Thus, the second 7-ray selection criterion was < w > < 1.2. This 
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set of cuts was found to be optimal for spectrum studies (Aharonian et al. 1999b, 1999c). 
These loose analysis cuts provide a Crab Nebula 7-ray rate of 83 7s/lir at SZA (less that 
25°) for the 5-IACT system. The corresponding energy threshold of the 7-rays is about 
500 GeV. For the LZA data the looser orientation cut of < 0.1 [deg^] was used because 
of the lower accuracy of the arrival direction reconstruction for the 7-ray showers. In order 
to improve the cosmic ray rejection in observations at LZA an additional parameter, mean 
scaled Length, < I > was used, which is defined by analogy with < w > (Konopelko et al. 
1999b). These two parameters, < w > and < I >, can be used for calculating a Mahalanobis 
distance, MD (Mahalanobis 1963), in two-dimensional space as 

MD = ((1- < w >r/al^^ + (1- < i >)V<i>)'/' (2) 

where cr<^> and c^iy are the standard deviations for the corresponding distributions of 
< w > and < / >. The optimum value of the MD cut for LZA is found to be 1.5. Note 
that this analysis improves the enhancement factor by ~ 30% (it gives ~ 50% acceptance 
of 7-rays) in observations at LZA, whereas it gives only marginal improvement for the data 
taken at SZA. The Crab Nebula 7-ray rate in observations at LZA (60°) is about 16 7's/hr 
with a corresponding energy threshold of ~5 TeV. Note that SZA observations give a 7-ray 
rate at high energies (above 3 TeV) of ~ 87s/hr. A summary of the data is shown in 
Table 2. 

The observations of the Crab Nebula have been made during 6 periods which differ 
in the system configuration, refiectivity of the mirrors, light refiection by the pixel funnels 
and camera protecting plate etc. All that afi^ects the hardware event rate of the telescope 
system, Rexp- These changes of system performance were implemented in the Monte Carlo 
simulations. Assuming the standard chemical composition of the primary cosmic rays 
(Wiebel, 1994) the calculated detection rates, Rmc, were adjusted to the measured rates 
(see Table 3). 
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The collection function of energy and zenith angle, for 7-ray showers have 

been inferred from Monte Carlo simulations (Konopelko et al. 1999a). The rms error of 
the energy determination is AE/E ~ 0.18. The Monte Carlo studies show that for a good 
energy resolution of 18% this approach does not distort the initial spectrum shape. The 
collection area for 7-rays rises very quickly in the energy range near the energy threshold 
of the telescope system, which is 500 GeV, whereas it is almost constant at the energies 
above > 3 TeV. Even slight variations of the trigger threshold could lead to noticeable 
systematic changes in the predicted spectral behavior in the energy range of ~ 0.5 — 1 
TeV. This effect leads to a noticeable probability for "sub-threshold" triggers. In addition, 
the trigger level for different camera pixels is slightly different even after very accurate 
adjustment of the high voltage using the cahbration laser runs. Measurements of the trigger 
setting for a number of camera pixels revealed variation in the trigger threshold of order 
10%. These variations were implemented into the simulations in order to estimate the 
corresponding systematic error of the energy spectrum at energies below 1 TeV. The fine 
tuning of the Monte Carlo simulations with respect to the lACT system data provided 
measurements of the flux of the cosmic ray protons in the energy range from 1.3 to 10 TeV 
(Aharonian et al. 1999d). The proton fluxes as measured by the HEGRA lACT system are 
perfectly consistent with the results of a bulk of satellite experiments held in this energy 
range. 

The procedure for the evaluation of the energy spectrum using the stereoscopic 
observations was discussed in detail by Aharonian et al. (1997); Hofmann (1997), and, 
more recently in Aharonian et al. (1999c). In the stereoscopic observations the impact 
distance of the shower axis to a system telescope can be measured with an accuracy < 10 
m. The energy E of a 7-ray shower is defined by interpolation over the "size" parameter S 
(total number of photoelectrons in Cerenkov light image) at a fixed impact distance R, as 
E = fMc(S, R, 6'), where 9 is the zenith angle and Imc is a function obtained from Monte 
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Carlo simulations. Note that the Monte Carlo simulations used here include the samphng 
of detector response in great detail (Hemberger 1998). The energy distribution for the ON- 
and OFF-source events, after the orientation and shape image cuts, were histogrammed 
over the energy range from 500 GeV to 30 TeV with 8 bins per decade. The 7-ray energy 
spectrum was obtained by subtracting ON- and OFF-histograms and dividing the resulting 
energy distribution by the corresponding collection area and the 7-ray acceptance. In the 
present Crab Nebula analysis the energy spectrum measurements were extended up to large 
zenith angles (65°). The data were processed independently for each of the four zenith 
angle bins: (0° - 25°), (25° - 40°), (40° - 50°), (50° - 65°). The corresponding effective 
collection areas as well as the cut efficiencies were calculated as a function of the zenith 
angle. First, the energy spectra were derived for all zenith angle bins independently. Note 
that the spectra evaluated at different zenith angles are in a good agreement. For the final 
energy spectrum the different zenith angle bins were joined according to the prescription: 

4 

dJl^/dE = J2w,{dj;/dE)ME' - El,), 
3=1 

Wj ^tj/to, i ^ l,n; (3) 

where dJl^/dE, {dJl^)j/dE are the differential energy spectra at energy as measured over 
all zenith angle ranges, and for the particular zenith angle bin (j), respectively, i?^^ is an 
estimated energy threshold for the zenith angle bin j, tj is the observation time for the j-bin 
on the zenith angle, and Iq is the total observation time. The first three zenith angle bins 
were joined using the time dependent weights wj = tj/to. Finally the spectrum measured 
in the zenith angle range of (0° — 50°) was combined with the spectrum derived from large 
zenith angle data, (50° — 65°) using the weights based on the estimated statistical errors 
for both spectra. Such procedure takes into account the advantageous 7-ray rate in LZA 
observations at high energies. 

The statistics of the 7-rays from the Crab Nebula provides a measurement of the 
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energy spectrum up to a few tens of TeV. However, detection of Cerenkov light images 
with extremely large amplitudes - several thousands of ph.e. - is complicated by the 
nonlinearity in the PMT response as well as by the saturation in the 8 bit Flash- ADC 
readout. Measurements of the photomultiplier response under high light loads over the 
extended sample of the EMI 9073 PMTs gave a cahbration function which was used to 
correct the image amplitudes. The readout of the HEGRA lACT is based on the sampling 
of Cerenkov hght time impulse by the 16 FADC bins of ~ 8 ns each (Hess et al. 1998). 
The time pulses from the air showers with a full width at half maximum of a few ns were 
widened using an electronic scheme in order to fit into several FADC bins for the accurate 
measurement of the time profile. The smoothing of the FADC signal was unfolded back to 
the impulse, which almost always fits 2 FADC bins. The calibrated amphtude, summed 
over two FADC bins, is used as a measure of the pixel signal. For the high energy air 
showers the FADC signals run into saturation and the simple unfolding procedure fails. For 
such pulses the initial amplitude is reconstructed using the additional calibration function 
obtained by simultaneous measurements of light flashes with FADCs and a 14 bit ADC. 
This procedure drastically extends the dynamic range of the FADC readout. 

To avoid the saturation problem one might only use images detected from air showers 
at large impact distances from the telescope system (e.g. beyond 150 m). The size of 
these images is very small even for high energy events because of the low Cerenkov light 
density far off the shower axis. However these images are very often truncated by the 
camera edge and do not allow a proper reconstruction of the shower impact point and of 
the shower energy. This effect becomes less important in observations at LZA because of 
the high shower maximum height (the images shrink to the camera center). In the present 
analysis the maximum impact distance of the shower core from the center of the system was 
extended up to 400 m for observations at LZA. Observations at LZAs permit measurements 
of the energy spectrum far beyond 10 TeV. The images of 7-ray air showers observed at 



-13- 



LZAs have small Size and are not influenced by the saturation effect. 

4. Results 

We have observed the Crab Nebula extensively in two observational seasons with the 
HEGRA lACT system. The HEGRA system of 5 lACTs currently has a sensitivity which 
allows the detection of a ^5cr signal from the Crab Nebula within 1 hr of observation 
time (see Figure 1). The integral 7-ray fluxes measured during the different observational 
periods are consistent within the estimated statistical and systematic error (see Table 3). 
The differential energy spectrum of the Crab Nebula has been derived from the HEGRA 
data for two observational campaigns using recently developed advanced techniques for 
the measurements of the spectrum using stereoscopic data taken at small and large zenith 
angles. The Crab Nebula differential energy spectrum derived from SZA data matches 
quite well the spectrum derived at LZA (see Figure 3) . The 7-ray rate measured at energies 
above 10 TeV in observations at LZAs exceeds the corresponding rate measured at SZA 
by a factor of 3. The LZA data are not affected by saturation effects. At 3.7(j confidence 
level 27 7-ray events from the Crab Nebula were detected in the highest energy bin from 
17.8 to 23.7 TeV (see Figure 4). One may expect that a number of 7-ray events in the 
highest energy bin are spilled over from the lower energies. However, given the good energy 
resolution of 20% and the power law energy spectrum, such effect is very small and is 
compensated almost by the backwards influx of the 7-rays from the energies above (see 
e.g., Aharonian et al. 1995). Finally, the simulations show that spilling over of low energy 
7-rays does not influence the resulting fluxes measured at the upper end of the power law 
spectrum. 

The analysis for the different system conflgurations as well as for different trigger 
threshold values gives a differential energy spectrum of the Crab Nebula measured at zenith 
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angles up to 60° 

dJ^/dE = (2.79 ± 0.02 ± 0.5) ■ 10-^(p^)-2-59±° °3±o °5phm-2s-iTeV-i (4) 

The statistical and systematic errors are also given. The final Crab Nebula spectrum as 
measured by the HEGRA collaborations is shown in Figure 4. The measured 7-ray fluxes 
are given in Table 4. The Crab Nebula energy spectrum is best fitted by a pure power law 
in the energy range 1-20 TeV. It does not exclude a possible slight steepening of the energy 
spectral usually predicted by inverse Compton modeling of TeV 7-ray emission. A fit with 
a logarithmic steepening of the power law spectrum gives the following result 

dJ^/dE = (2.67 ± 0.01 ± 0.5) ■ io-7(_|_)-2.47±o.i±o.o5-(o.ii±o.io) iog(E)pj^ ^-2 g-i Tev-i (5) 

Such a fit indicates the slight fiattening of the spectrum at low energies as predicted by the 
IC calculations. However, the change of the energy spectrum slope is within the current 
statistical and systematic errors, and the data for the overall energy range are consistent 
with a simple power law fit in 0.5-20 TeV. The HEGRA Crab Nebula data match well the 
recent 20 TeV data published by the CANGAROO group (Tanimori et al. 1998), and are 
consistent with a flat power law index of ~ 2.5 beyond 20 TeV. The compilation of the 
world data is given in Figure 5. All data are consistent within statistical and systematic 
errors, except possibly for the Tibet data which show relatively higher fluxes. 

5. Astrophysics implications 

The TeV energy spectrum of the Crab Nebula as measured by the HEGRA system of 
lACTs is consistent with the expectations for the TeV 7-ray emission from pulsar-driven 
Nebulae (plerions). According to this scenario the ultra relativistic electrons, accelerated 
in the pulsar wind shock, produce TeV 7-rays through the IC scattering with soft photons 
within the Nebula. The predicted fluxes of TeV 7-rays rely on the spatial distribution of 
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the magnetic field within the Nebula as determined by the parameter a (ratio of energy 
density of magnetic field to the particle energy density) and/or by the average magnetic 
field < B > in the optical nebula. The HEGRA data are shown in Figure 6 together with 
predicted spectra using two SSC models of TeV 7-ray emission (de Jager et al. 1996; Atoyan 
& Aharonian 1996). One may conclude that both models fit the HEGRA data rather well. 

According to the calculations of the TeV 7-ray emission by de Jager et al. (1996) the 
7-ray flux from the Crab Nebula at TeV energies constrains the choice of the parameter a. 
The IC spectrum computed by de Jager et al. (1996), assuming for the parameter a a value 
~ 0.003 gives a good fit to the HEGRA data. This value of the parameter a corresponds 
to the best-fitting magneto-hydrodynamic (MHD) solution of electron propagation in the 
Crab Nebula as found by Kennel & Coroniti (1984). 

In another approach, assuming the spatial distribution of magnetic field in the Crab 
Nebula, one can determine the average magnetic field < B > in the optical nebula (Gould 
1965). According to the calculations of Atoyan & Aharonian (1996), made within the 
framework of the MHD model of Kennel & Coroniti (1984), the average magnetic field < B > 
is determined by the TeV 7-ray flux as < B >oc J°-^ • 10~^G where J^, phcm~^s~^TeV~^ 
is a differential flux at 1 TcV. Thus the HEGRA spectrum gives an average magnetic field 
strength < B >~ (1.7 ± 0.3) ■ 10~'^G. This value is consistent with the estimate derived by 
Hillas et al. (1998) from the Crab Nebula date taken with the 10 m Whipple telescope. 

The IC energy spectrum of 7 rays from the Crab Nebula, measured in the energy 
range from 1 TeV to 10 TeV, is hkely to be a power law dJ^/dE oc E"" = E-^-^ At the 
same time, for the energy range above 10 TeV calculations predict gradual steepening with 
a ~2.7 and 2.9 at 10 and 30 TeV, respectively. That is due to both the energy loss of the 
ultra high energy electrons by fast synchrotron cooling and the Klein-Nishina effect in the 
cross-section of the inverse Compton scattering. Atoyan & Aharonian (1996) and Bednarek 
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& Protheroe (1997) have shown that 7r°-decay 7-ray fluxes, due to the relativistic protons 
accelerated in the Crab Nebula, may noticeably contribute at energies above 10 TeV. 
However the HEGRA Crab Nebula data expanded up to 20 TeV are still consistent with the 
pure IC spectrum. To assess the contribution of 7r°-produced 7-rays from the Crab Nebula 
measurements above 30 TeV are needed. Note that the LZA technique could help to 
perform such observations. 

The predicted IC 7-ray spectrum of the Crab Nebula is rather fiat in the energy range 
below ~ 1 TeV. It could be well approximated by dJ^/dE oc E~^ ° at 100 GeV. Detection of 
a gradual flattening in this energy range will prove the SSC scenario of TeV 7-ray emission. 
However the low energy points at the HEGRA Crab Nebula spectra (E-^ < 1 TeV) are 
strongly affected by possible systematic errors (~ 50%) and do not allow such a conclusion. 
Future observations of the Crab Nebula with the forthcoming low threshold high sensitivity 
Cerenkov detectors (see Catanese & Weekes 1999) will offer precise measurements in these 
energy range. 
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Fig. 1. — Distribution over Q"^ {9 is an angular distance of the true source position to the 
reconstructed source position) of ON (sohd hne) and OFF (dashed hne) events in 12 hrs of 
observations of the Crab Nebula with the HEGRA lACT system at zenith angles < 20°. 

Fig. 2. — Filled contour plot (gray scale) of the smoothed distribution of the reconstructed 
core shower positions in the sky (Dec. vs RA) in a 1" x 1° region centered around the position 
of the Crab. Superimposed are contours of the significance (the contours increment by one 
a). Number of total events: 27800; number of events in the peak: 180. All events have 
< w >< 1.2. 

Fig. 3. — The Crab Nebula energy spectrum as measured at small and large zenith angles 
using the HEGRA lACT system. The filled circles are for the observations at zenith angles 
up to 50°, the open circles are for the LZA data (60°). 

Fig. 4. — The Crab Nebula energy spectrum derived from the data taken over the zenith 
angle range from 0° to 65°. The power law fit of eqn. (4) is shown by the solid hne. The 
power law fit with logarithmically energy dependent slope of eqn. (5) is shown by the dashed 
curve. 

Fig. 5. — Differential spectrum of the Crab Nebula as measured by the HEGRA sys- 
tem of imaging air Cerenkov telescopes in comparison with the results of Whipple 
(Hillas et al. 1998), CANGAROO (Tanimori et al. 1998), CAT (Barrau 1998), and Ti- 
bet (Amenomori et al. 1999). The dashed curves show the upper and lower limits of the 
systematic error estimated for the HEGRA data. 

Fig. 6. — Spectral energy distribution of TeV 7-ray emission from Crab Nebula. The HEGRA 
data are shown by dots. Results of calculations using SSC model of emission made by Atoyan 
& Aharonian (1996) and de Jager et al. (1996) are shown by sohd and dotted-dashed curve, 
respectively. 
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Table 1. Summary of exposure times. 



N 


Period 


Number of tel-s 


Z.A. 


Time [hrs] 


1 


1997 Sep - 1998 Jan 


4 


0° - 50° 


45.6 


2 


1997 Oct - Nov 


3 


0° - 50° 


14.3 


3 


1998 Feb - March 


4 


0° - 50° 


22.6 




Total: 






82.5 


4 


1998 Sep - Oct 


4 


0° - 50° 


26.2 








50° - 65° 


5.4 


5 


1998 Oct - Nov 


5 


50° - 65° 


29.8 








50° - 65° 


10.4 


6 


1998 Dec 


5 


0° - 50° 


20.1 








50° - 65° 


8.1 




Total: 




0° - 50° 


76.1 








50° - 65° 


24.0 
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Table 2. Summary of the data. No. of 7-rays, cosmic rays (CRs), and the S/N ratio, in 

the standard deviations {a). 





Z.A. 


0° - 25° 


25° - 


40° 40° - 50° 


55° - 65° 




Time, hrs 


7.30 


11.34 


7.56 


5.44 




Raw 


592 


972 


327 


143 


No. of 7's 


L. cut 


577 


670 


244 


124 




MD cut 


450 


499 


157 


101 




Raw 


2964 


4495 


2772 


3353 


No. of CRs 


L. cut 


536 


846 


567 


658 




MD cut 


191 


268 


151 


114 




Raw 


7.3 


9.7 


4.3 


1.8 


S/N [a] 


L. cut 


14.2 


13.8 


6.6 


3.3 




MD cut 


15.6 


15.5 


7.3 


5.6 



Note. — Data for the system of 4 lACTs were analyzed using a ori- 
entation cut of < 0.05 deg^, and a loose mean scaled Width cut of 
< w >< 1.2 (L. cut), a Mahalanobis distance cut (MD< 1.5) as well as a 
cut without image-shape (Raw). 
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Table 3. The integral 7-ray flux above 1 TeV, ^^{> 1 TeV), in units of photon cm^^s"-'^, 
measured during the different observation periods. Rcxp and Rmc ai'e the measured and 
calculated hardware detection rates, respectively. The data correspond to a zenith angle 

range from to 50 degrees. 



Period 


Rexp 


R-MC 








1 


12.6 


14.0 


(1.68 ±0.04) 


• 10- 


-11 


2 


10.3 


11.2 


(1.71 ±0.05) ' 


■ 10- 


-11 


3 


10.7 


10.0 


(1.68 ±0.05) ■ 


■ 10- 


-11 


4 


11.5 


11.4 


(1.77 ±0.04) ■ 


• 10- 


-11 


5 


10.1 


10.0 


(1.84 ±0.04) 


■ 10- 


-11 


6 


11.8 


11.4 


(1.64 ±0.05) ■ 


■ 10- 


-11 



Note. — These results were obtained using 
the power law fitting of the differential 7-ray 
fluxes measured in the energy range from 1 to 
10 TeV for each observational period. 
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Table 4. The energy spectrum of the Crab Nebula, dJ-y/dE, in unites of photon cm~^s~ 
The statistical ((Jstat) and systematic ((Jstat) errors shown in the table are also given in 

photon cm~^s~^. 



E, TeV 


dJ^/dE 


Cstat 










7.44-10- 


■11 


5 25-10" 


-12 


(+2 31-1 79"! 10" 


-11 


0.87 


3.74-10- 


-11 


2.50-10" 


-12 


(+5.98-5.61)10- 


-12 


1.16 


1.95-10" 


■11 


1.39-10- 


-12 


(+1.94-1.94)10- 


-12 


1.54 


9.25-10- 


■12 


1.58-10" 


-13 


(+8.32-7.40)10- 


-13 


2.05 


4.25-10- 


-12 


1.03-10" 


-13 


(+3.40-2.97)10- 


-13 


2.74 


2.20-10- 


■12 


7.72-10" 


-14 


(±1.54)10-^3 




3.65 


9.78-10- 


13 


5.17-10- 


-14 


(±5.87)10-1^ 




4.87 


4.43-10- 


13 


3.25-10" 


-14 


(±2.66)10-^^ 




6.49 


2.32-10- 


13 


3.40-10" 


-14 


(±1.39)10-^^ 




8.66 


1.20-10- 


-13 


2.45-10" 


-14 


(±7.19)10-15 




11.55 


5.64-10- 


■14 


1.68-10- 


-14 


(±3.38)10-15 




15.40 


2.28-10- 


■14 


1.07-10" 


-14 


(±1.36)10-15 




20.54 


1.14-10- 


■14 


7.55-10" 


-15 


(±6.85)10-16 
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5.54 5.55 5.56 5.57 5.58 5.59 5.6 5.61 

RA (hours) 
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^ 10 

CO 

CM 

E 
o 

> 

CD 



LU 
T3 



LU 



-10 



10 



-11 



10 



-12 



10 




★ OO Whipple (1998) 

A CANGAROO (1998) 

CAT (1998) 
□ Tibet (1999) 

• HEGRA GT-System 



□ 



10 



A 



10 
E [TeV] 
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